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SUMMARY 


The  procedure  for  designing  a  detour  phase  computer  generated 
holograms  is  described  with  particular  reference  to  a  scale  and 
rotation  invariant  transformation  filter  known  as  a  Mellin  transform 
filter.  Experimental  results  are  presented  which  demonstrate  the 
main  features  of  the  devices. 
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1  INTRODUCTION 


Coherent  optical  systems  possess  the  great  qualities  of  high  data  through¬ 
put  with  parallel  processing.  In  the  last  few  years  several  areas  which  utilise 
these  qualities  have  developed  and  are  becoming  increasingly  important  with  the 
development  of  high  resolution  holographic  film  and  real  time  coherent  light 
valves.  Phase  masks  which  were  not  possible  to  achieve  using  purely  optical 
systems  are  now  being  developed  using  the  various  techniques  of  computer 
generated  holography^^”^®^.  Likewise  real  time  optical  correlator  systems 
have  been  demonstrated  using  the  new  generation  of  optical  light  valves.  The 
two  main  optical  correlators  that  are  usually  demonstrated  are  the  frequency 
plane  correlator  and  the  Joint  transform  correlator.  Both  of  these  systems 
exhibit  spatial  invariance  but  are  sensitive  to  scale  and  rotation.  In  1983 
Saito  et  al^^)  described  a  scale  and  rotation  invariant  real  time  optical 
correlator  which  employed  a  Mellin  transform  phase  mask.  The  purpose  of  this 
memorandum  is  to  describe  the  methods  employed  to  create  this  Mellin  transform 
filter  and  to  discuss  some  of  its  limitations.  During  this  exercise  it  is 
hoped  that  the  reader  will  obtain  some  insight  into  the  field  of  computer 
generated  holograms. 


2  CONTINUOUS  PHASE  DESCRIPTION  FOR  A  MELLIN  TRANSFORM  FILTER;  THEORY 


Consider  the  Fourier  transform  properties  of  a  simple  lens.  The  light 
intensity  in  the  front  focal  plane  of  a  lens,  F(u,v),  is  related  to  the  light 
intensity  in  the  rear  focal  plane,  f(x,y),  by  the  relationship 


F(u,v) 


•Jh 


y)  exp(-jkxu/f)  exp(-jkyv/f)  dxdy 


(1) 


where  f  is  the  focal  length  of  the  lens  and  k  •  A  is  the  wavelength  of 

the  light.  In  our  particular  case  we  consider  an  Initial  input  function  which 
has  only  amplitude  information,  the  phase  content  of  the  function  f(x,y)  is 
supplied  by  a  phase  plate  that  is  coplanar  with  this  input  function,  thus 
f(x.y)  ■  gCx.y)  exp(J<f(x,y))  and 


F(u,v) 


fh 


y)  exp( j(t(x.y)  ~  kxu/f  -  kyv/f))  dxdy 


(2) 


By  applying  the  method  of  stationary  phase  to  evaluate  this  integral  it  can  be 
shown  that  the  point  (x,y)  of  plane  f(x,y)  is  mapped  into  the  point  (u,v)  of 
plane  F(u,v)  by  the  following  relationships 


Xf  3$  , 


2lI  li  (x  v) 
2n  3y 


(3) 


For  a  continuous  phase  computer  generated  hologram  we  require  that  the  func¬ 
tions  u  and  V  have  continuous  partial  derivatives  in  a  simply  connected  region 
of  plane  f(x,y).  This  Implies 


3y  3x  3y  3x 


(A) 


In  most  cases  of  computer  generated  holograms  the  mapping  that  is  required 
is  defined  and  it  is  necessary  to  calculate  the  appropriate  phase  function 


^(x»y)  using  equation  (3)  (subject  to  the  continuous  phase  conditions  of 
equation  (4)).  In  Appendix  A  the  case  Is  considered  when  equation  (4)  Is  not 
satisfied. 


In  the  present  case  we  choose  the  laapplng 

0/u\  /£n(r)*s\  /£n(x^+y^)  •  s/2\ 

\v/  \  -0s  /  \-tan”'(y/x)  •  s/ 


(5) 


Careful  evaluation  of  equation  (4)  confirm  that  this  function  does  Indeed 
produce  a  continuous  phase  function  and  solving  the  two  differential  equations 
(3)  yield  the  final  form  of  the  phase  function  (Ref  1) 
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where  s  is  a  scaling  factor  and  H(-x)  Is  the  unit  step  function  which  has  the 
following  definition 


H(-x)  -  1  X  <  0 

H(-x)  -  0  X  >  0 

THE  BINARY  FRAUNHOFER  HOLOGRAM  AND  THE  MELLIN  TRANSFORM  FILTER 


(7) 


A  binary  hologram  consists  of  transparent  dots  on  an  opaque  background: 
the  transmittance  of  the  medium  being  either  one  or  zero  at  all  places  on  the 
transparency^^ In  order  that  a  complex  amplitude  f(x,y)  can  be  realised 
on  the  hologram  we  must  consider  how  phase  information  can  be  recorded.  There 
are  three  ways  in  which  the  phase  of  a  light  wave  can  be  Influenced:  retarda¬ 
tion  while  travelling  through  a  dielectric,  phase  jump  at  reflection  and  detour 
phase.  The  mechanism  of  detour  phase  can  be  visualised  using  the  description 
of  grating  diffraction.  In  the  first  order  diffraction  pattern  light  from 
adjacent  slits  have  a  relative  2x  phase  shift,  this  Is  the  detour  phase  and 
If  the  relative  centres  of  the  slits  are  not  correctly  placed  It  is  clear  that 
an  alternative  value  for  the  detour  phase  can  be  realised. 

In  the  synthesis  of  the  complex  filter  function  it  Is  assumed  that  the 
cells  are  of  such  a  size,  (Ax)^,  that  the  function  f(x,y)  will  almost  be 
constant  across  the  cell 


f(x,y)  =s  f(nAx,  mAy)  -  fnm  -  Anm  exp(j(|>nm) 

The  cells  may  also  be  so  small  that  the  envelope  function  that  they  add  to 
the  waveform  Is  unimportant.  Several  ways  of  implementing  the  phase  and 
amplitude  information  have  been  tried.  For  the  Mellin  transform  filters  In 
this  memorandum  the  schemes  shows  In  Figure  1  have  been  used.  The  simplest 
case  consists  of  a  slit  whose  width  Is  half  of  the  cell  size  and  whose  height 

Is  normalised  so  that  the  maximum  amplitude  required  by  the  function  corres¬ 
ponds  exactly  to  the  cell  height.  The  position  of  this  element  relative  to 

the  centre  of  the  cell  Is  a  measure  of  the  detour  phase  required  of  this 

element.  One  disadvantage  of  this  scheme  Is  that  there  is  a  slight  possibility 
of  the  contribution  from  adjacent  cells  overlapping.  This  may  be  overcome  by 
folding  the  aperture  round  within  the  cells  so  that  we  effectively  either  add 
or  subtract  2x  from  the  detour  phase  of  the  contribution  outside  the  cell  to 
bring  It  back  within  the  cell  walls  (see  Figure  1). 


Il 


When  diffracted  light  from  a  regular  array  of  cells  of  dimension  dx  Is 
collected  by  a  lens  of  focal  length  f,  the  diffraction  pattern  seen  In  the 
focal  plane  of  this  lens  Is  a  regular  square  array  of  points  where  the 
separation  of  the  points  Au  Is  given  by 

Au  •  f  =  Av  (9) 


We  have  chosen  to  map  the  (x,y)  plane  to  the  (u,v)  plane  In  such  a  way  that  a 

rotation  of  an  object  about  the  origin  by  an  angle  6  in  the  (x,y}  plane  will 

correspond  to  a  linear  shift  of  Os  In  the  v  direction  of  the  (u,v)  plane. 

We  choose,  therefore,  the  size  of  the  cells  In  the  hologram  to  be  such  that 

the  separation  of  the  diffraction  orders  Av  In  the  (u,v)  plane  to  correspond 

to  an  Integer  multiple  of  2vs,  the  linear  shift  introduced  In  the  v  direction 
due  to  complete  rotation  of  the  Image  in  the  (x,y)  plane. 

In  this  way  the  diffraction  order  spots  will  give  a  direct  calibration 
of  6  and  £n(r)  In  the  (u,v)  field.  Thus 


Ay 


Xf  Xf  Xf 

■■■  ■  s  —  s  _ 

27rns  Au  Av 


(10) 


all  linear  dimensions  are  measured  In  the  same  units.  Next  we  choose  the 
separation  of  the  orders  in  the  (u,v)  plane  to  be  5  mm,  a  convenient  size  to 
image  on  our  TV  system.  Thus  Av  ■  5  ■  2’rs,  The  parameter  s  therefore 
takes  the  values  (5/2r)  and  (1/2*)  respectively  for  the  two  conditions 
(n»l)  and  (n«5)  considered  In  the  experiments  (using  equation  (10)),  For 
X  •  514.5  X  10"^  mm  and  f  ■  1000  mm  the  above  choice  of  Av  corresponds  to  a 
cell  size  In  the  hologram  of  102.9  pm  square. 

The  Mellln  transform  filter  that  would  be  produced  using  these  parameters 
was  simulated  using  the  'basic*  program  Included  In  Appendix  B.  It  is  clear 
from  these  programs  that  an  amplitude  term  has  been  Included  In  the  plot.  The 
Mellln  filter  Is  a  phase-only  mask  and  the  amplitude  term  Is  Included  so  that 
a  linear  magnification  of  the  Input  image  will  result  in  a  corresponding  shift 
In  the  output  pattern  but  will  have  no  resulting  change  of  Intensity  of  the 
transform  Image.  In  the  absence  of  such  a  term  a  magnified  Input  Image  would 
have  a  transformed  Image  that  was  much  brighter.  This  is  not  desirable  when 
the  detecting  system  may  have  a  limited  dynamic  range  and  where  one  would  wish 
to  define  a  suitable  operating  point  that  Is  Invariant  with  respect  to  scale. 
As  an  Image  intensity  will  scale  according  to  an  r^  term  with  respect  to 
Iniform  illumination,  the  amplitude  term  must  vary  as  r"^  (the  amplitude  in 
our  plots  Is  set  arbitrarily  to  unity  for  r  <  1  mm,  and  scaled  as  1/r  for 
r  >  1  mm).  The  resulting  filters  are  shown  In  Figure  2.  The  final  holograms 
were  produced  using  the  Electromask  facility  at  RSRE.  This  facility  has  a 
minimum  slit  of  2.6  pm  amd  a  maximum  slit  of  1.5  mm.  The  positioning  accuracy 
of  the  slit  Is  0.25  pm  over  a  60  nun  field.  The  size  of  transform  filter  we 
can  reasonably  produce  Is  limited  by  the  writing  speed  of  the  electromask 
(4  k/hr  for  chrome  masks,  16  k/hr  for  emulsion  masks).  For  our  particular 
device  the  write  time  of  between  2  and  8  hours  was  required  for  the  3  x  10^ 
Images. 
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EXPERIMENTAL  RESULTS 


The  first  Hellin  transform  filter  (MELLIN  No  1)  has  been  used  to  demon¬ 
strate  the  in(r)-6  transformation  in  various  situations.  All  of  the  experimental 
results  were  taken  using  the  system  depicted  in  Figure  3.  Light  from  an  Argon 
ion  laser  is  spatially  filtered  and  collimated  on  to  the  transform  filter  which 
is  in  the  near  focal  plane  of  a  1  m  focal  length  achromatic  doublet  lens.  The 
lens  collects  the  light  and  focusses  it  on  to  the  active  surface  of  a  lowlight 
TV  camera  which  lies  in  the  front  focal  plane  of  the  lens.  The  image  is  then 
displayed  on  a  suitable  monitor.  Amplitude  modulated  input  data  is  placed  just 
behind  the  Mellin  filter  in  the  collimated  light. 

The  active  area  of  the  filter  has  a  10.29  mm  radius.  Illumination  of  the 
whole  of  this  area  corresponds  to  a  band  in  the  in(r)-6  plane  as  depicted  in 
Figure  4.  This  defines  the  region  of  the  transform  plane  (u,v)  that  meaningful 
results  can  be  measured.  The  upper  edge  of  the  Illuminated  region  of  the  (u,v) 
plane  corresponds  to  r  ■  10.29  (in(r)  ■  2.33  •  0.74it).  Correspondence  between 
the  theory  and  experiment  can  be  seen  to  be  very  good. 

When  the  aperture  has  a  restricted  outer  radius  the  upper  edge  of  the 
image  in  the  transform  plane  lowers  (see  lower  part  of  Figure  5).  If  the  lower 
radius  of  the  image  is  also  restricted  so  that  the  filter  is  illuminated  by  a 
circular  annulus  (Figure  6)  then  the  transformed  image  becomes  a  narrow  band 
in  the  (u,v)  plane  whose  width  is  given  by  in(b/a)  where  a  and  b  are  the  inner 
and  outer  radii  respectively.  If  this  circular  annulus  is  further  restricted 
so  that  only  a  half  of  the  aperture  is  transparent  (Figure  7)  the  image  becomes 
restricted  in  the  6-dlrection  also.  Figure  7  demonstrates  the  effect  of  rota¬ 
tion  of  the  circular  half  annulus  on  the  transformed  image.  It  is  apparent 
that  the  major  effect  is  not  in  the  fine  detail  of  the  image  but  is  merely  a 
translation  of  the  same  image  along  the  6  axis.  This  overcomes  one  of  the 
problems  of  the  space  invariant  transform  -  namely  that  of  Invariance  to 
rotation.  This  effect  is  further  demonstrated  in  Figure  8  with  a  full  square 
annulus.  In  this  case  the  peak  features,  corresponding  to  the  corners  of  the 
square,  shift  by  ir/4  in  the  transform  plane  upon  a  rotation  of  the  illuminating 
image  by  a  similar  amount.  Note  also  in  this  figure  the  multiple  diffraction 
pattern  associated  with  the  cell  diffraction  spots.  This  pattern  is  a  repli¬ 
cated  version  of  the  pattern  chat  one  would  observe  at  the  centre  of  the  trans¬ 
form  plane  if  a  single  cell  of  dimensions  4x^  is  placed  in  the  input  plane,  and 
is  useful  in  confirming  the  orientation  of  the  square  annulus  with  respect  to 
the  Mellin  axes. 

When  an  image  is  no  longer  on  axis  and  does  not  even  Include  the  origin 
in  the  x,y  plane,  Che  pattern  in  the  transform  plane  is  very  different.  For  a 
circular  aperture  (Figure  5)  the  pattern  is  an  oblaced  circle  whose  size  and 
distortion  depends  on  its  position  with  respect  to  the  (x,y)  origin.  If  a 
more  complex  shape,  eg  the  square  annulus,  is  taken  off  axis  its  pattern 
depends  strongly  on  the  orientation  of  the  square  with  respect  to  the  Mellin 
axes.  Two  such  cases  are  shown  in  Figure  9  and  Figure  10  for  two  orientations 
of  Che  square  which  differ  by  a  rotation  of  x/4.  The  patterns  are  obviously 
very  different  although  the  same  shape  is  employed  but  Che  correspondence 
between  experiment  and  theory  is  still  excellent.  This  demonstrates  the  fact 
that  Che  Mellin  transform  can  only  show  rotation  Invariance  with  respect  to  its 
own  centre  -  not  the  centre  of  the  concerned  object  that  illuminates  the  filter. 
This  is  a  particularly  Important  point  with  respect  to  use  of  the  filter  in  a 
correlator  system  -  obviously  targets  of  interest  must  be  aligned  before 
analysis  if  this  scheme  is  to  be  used. 


In  this  section  experimental  data  has  only  been  displayed  for  the  first 
Mellin  transform  filter  (Mellin  No  1).  An  Improved  version  of  this  filter 
was  constructed  (Mellin  No  3)  using  the  alternative  method  of  CGH  construction 
shown  in  Figure  1.  The  filter  Itself  looks  more  aesthetically  pleasing  to  the 
eye  (see  Figure  2)  showing  no  large  gaps  in  the  filter  and  displaying  a  more 
continuous  form*  Experimentally  however  there  is  only  marginal  Improvement 
in  performance  using  this  filter  and  so  no  data  is  presented.  This  lack  of 
significant  improvement  is  probably  due  to  the  sampling  conditions  of  the 
object  (ratio  of  the  cell  size  to  object  size)  remaining  the  same  for  the 
two  filters.  An  alternative  filter  (Mellin  No  2)  hoped  to  overcome  some  of 
the  problems  incurred  due  to  the  near  undersampling  of  the  phase  function  by 
the  previous  two  holograms.  The  phase  function  was  expanded  by  a  factor  of 
5  and  this  resulted  in  a  contraction  of  the  output  image  of  a  similar  factor 
(Figure  11).  In  Figure  1 1  we  note  that  the  6  field  does  not  appear  to  be 
symmetrical  about  the  diffraction  orders.  At  first  this  may  appear  to  be 
incorrect  but  if  the  form  of  the  phase  function  is  reconsidered  (equation  (6)) 
it  is  noted  that  Images  in  the  2nd  and  3rd  input  quadrants  (negative  x)  have 
an  additional  linear  phase  term.  This  linear  term  acts  like  a  single 
sideband  modulator  shifting  the  field  produced  by  these  two  quadrants  to  an 
asymmetric  position  about  the  final  diffraction  pattern.  The  magnitude  of 
this  shift  is  such  that  the  field  for  positive  x  and  negative  x  meet 
precisely.  In  practice  this  effect  is  of  little  consequence  as  it  only 
produces  a  linear  shift  of  the  output  plane. 

5  CONCLUSIONS 

The  use  of  a  detour  phase  computer  generated  hologram  to  provide  a  phase 
plate  which  can  produce  a  given  transformation  has  been  demonstrated.  The 
filter,  when  used  in  conjunction  with  the  appropriate  coherent  light  and  lens 
of  correct  focal  length,  has  been  shown  to  produce  an  output  image  whose  shape 
is  independent  of  magnification  and  rotation  about  the  transform  filter  origin. 
This  filter  may  be  of  significant  use  in  particular  areas  of  optical  signal 
processing  where  the  direction  of  the  image  is  known  in  advance  (otherwise  a 
raster  type  scan  must  be  employed).  Of  the  three  filters  most  of  the  experi¬ 
mental  data  displayed  in  this  memorandum  was  obtained  using  Mellin  No  1 .  An 
improved  version  of  the  filter  using  a  modified  plotting  routine  (Mellin  No  3, 
Figures  1  and  2)  showed  no  significant  advantage  over  the  first  type.  It  is 
true  that  the  improved  filter  is  aesthetically  better  with  regard  to  the  absence 
of  obvious  errors  in  the  phase  description  but  the  limited  number  of  cells  that 
describe  particular  images  still  results  in  some  distortion  of  the  images, 
especially  with  regard  to  the  0  and  x  positions  in  the  transform  plane. 
Similarly  the  modified  plotting  routine  that  described  the  filter  with  a  slower 
varying  phase  characteristic  shared  no  significant  Improvement.  This  once 
more  was  probably  due  to  a  similar  number  of  cells  still  being  used  to  describe 
the  same  images.  The  major  experimental  difference  between  Mellin  No  2  and 
Mellin  Nos  1,  3  was  the  reduction  in  the  output  plane  that  was  usefully  employed 
(Figure  11). 

We  have  chosen  to  Implement  these  filters  using  binary  computer  generated 
holographic  techniques.  As  only  phase  only  functions  are  used  in  these  filters 
it  may  be  more  appropriate  to  develop  future  filters  using  alternative  tech¬ 
niques,  eg  Referenceless  On-Axis  Complex  Hologram  (ROACH)'^\  These  techniques 
have  the  advantages  of  high  efficiency  and  that  of  the  transformation  being 
directed  on-axis  and  not  to  the  first  order  diffraction  pattern. 
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FIGURE  1  DEFINITION  OF  CELL  APERTURES  FOR  REALIZING  CO,.  - 

GENERATED  HOLOGRAMS  USING  THE  DETOUR  PHASE  .Vi.  .  G.. 
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EXPERIMENTAL  ARRANGEMENT  USED  TO  CONFIRM 
OPERATION  OF  MELLIN  TRANSFORM  FILTERS 
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FIGURE  4 

THE  MELLIN  TRANSFORM  FILTER  (MELLIN  #1)  FULLY  ILLUMINATED.  COMPARISON 
OF  EXPERIMENTAL  DATA  WITH  THEORY.  EACH  SQUARE  HAS  DIMENSIONS  OF  ir/3  IN 
EACH  DIRECTION,  fin(r)  IS  PLOTTED  VERTICALLY,  6  HORIZONTALLY.  THE  DOTS 
CORRESPOND  TO  THE  DIFFRACTION  PATTERN  FROM  THE  CELL  SIZE  USED  IN  THE 
COMPUTER  GENERATED  HOLOGRAM 


THEORY  tn(r),  0  THEORY  x.  y  EXPERIMENT  (.n{r), 


THE  MELLIN  TRANSFORM  FILTER  (MELLIN  ^1)  ILLUMINATED  CY  A  CiiU'UI 
APERTURE  OFF  AND  ON  AXIS.  COMPARISON  OF  EXPERIMENTAL  DATA  Wl Tl : 
THEORY.  EACH  SQUARE  HAS  DIMENSIONS  OF  r.'3  IN  EACH  DIRECTION.  in(w  l: 
PLOTTED  VERTICALLY,  (;  HORIZONTALLY.  THE  DOTS  CORRESPOND  TO  iHE 
DIFFRACTION  PATTERN  FROM  THE  CELL  SIZE  USED  IN  THE  COMPUTER  GENERATED 
HOLOGRAM 


THEORY  fin{r),fl  THEORY  x,  y  EXPERIMENT  Bn{r),  0 


THE  MELLIN  TRANSFORM  FILTER  (MELLIN  #1)  ILLUMINATED  BY  A  CIRCULAR 
ANNULUS  ON  AXIS.  COMPARISON  OF  EXPERIMENTAL  DATA  WITH  THEORY.  EACH 
SQUARE  HAS  DIMENSIONS  OF  v/S  IN  EACH  DIRECTION.  en(r)  IS  PLOTTED 
VERTICALLY,  6  HORIZONTALLY.  THE  DOTS  CORRESPOND  TO  THE  DIFFRACTION 
PATTERN  FROM  THE  CELL  SIZE  USED  IN  THE  COMPUTER  GENERATED  HOLOGRAM 


,  S'-,' 


THEORY  ln(r),0 


THEORY  X,  y 


EXPERIMENT  vii(; 


Mm¥ 


Ir  “  ^')w- 


«  tOOaV, 


FIGURE  7 
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THE  MELLIN  TRANSFORM  FILTER  (MEI.LIN  tl)  ILLUMINATED  DY  A  HALF  CIRCULAR 
ANNULUS  ON  AXIS  AT  VARIOUS  ORIENTATIONS.  COMPARISON  OF  EXPERIMENTAL 
DATA  WITH  THEORY.  EACH  SQUARE  HAS  DIMENSIONS  OT  .r'S  IN  EACH  DIRECTION 
fn(r)  IS  PLOTTED  VERTICALLY,  ('  HORIZONTALLY.  THE  DOTS  CORRESPOND  TO  THE 
DIFFRACTION  PATTERN  FROM  THE  CELL  SIZE  USED  IN  THE  COMPUTER  GENERATED 
HOLOGRAM 
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THE  MELLIN  TRANSFORM  FILTER  (MELLIN  #1)  ILLUMINATED  BY  A  SQUARE 
ANNULUS  OFF  AXIS.  COMPARISON  OF  EXPERIMENTAL  DATA  WITH  THEORY.  EACH 
SQUARE  HAS  DIMENSIONS  OF  7r/3  IN  EACH  DIRECTION.  (n(r)  IS  PLOTTED 
VERTICALLY,  0  HORIZONTALLY.  THE  POTS  CORRESPOND  TO  THE  DIFFRACTION 
PATTERN  FROM  THE  CELL  SIZE  USED  IN  THE  COMPUTER  GENERATED  HOLO(iRAM 
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FIGURE  10 


THE  MELLIN  TRANSFORM  FILTER  (MELLIN  -'-‘1)  ILLUMINATED  BY  A  SOU'NRE 
ANNULUS  OFF  AXIS.  COMPARISON  OF  EXPERIMENTAL  DATA  WITH  THEORY.  EACH 
SQUARE  HAS  DIMENSIONS  OF  7t/3  IN  EACH  DIRECTION.  !n(r)  IS  PLOTTED 
VERTICALLY,  0  HORIZONTALLY.  THE  DOTS  CORRESPOND  TO  THE  DIFFRACTION 
PATTERN  FROM  THE  CELL  SIZE  USED  IN  THE  COMPUTER  GENERATED  HOLOGRAM 


THE  MELLIN  TRANSFORM  FILTER  (MELLIN  #2)  FULLY  ILLUMINATED.  COMPARISON 
OF  EXPERIMENTAL  DATA  WITH  THEORY.  EACH  SQUARE  HAS  DIMENSIONS  OF  5ir/3  IN 
EACH  DIRECTION.  en(r)  IS  PLOTTED  VERTICALLY,  0  HORIZONTALLY.  THE  DOTS 
CORRESPOND  TO  THE  DIFFRACTION  PATTERN  FROM  THE  CELL  SIZE  USED  IN  THE 
COMPUTER  GENERATED  HOLOGRAM 


APPENDIX  A 


SOME  ALTERNATIVE  GEOMETRIC  TRANSFORMATIONS  USING  COMPUTER  GENERATED  HOLOGRA.MS 


In  the  text  we  have  discussed  the  stationary  phase  method  of  realising 
a  geometric  transform  where  the  phase  function  has  a  continuous  form.  In 
general  this  may  not  be  always  true  and  we  must  consider  how  a  hologram  may 
be  constructed  when  this  is  not  the  case.  One  solution  is  to  subdivide 
the  hologram  into  8ubholograras(2,10);  where  the  phase  variation  within  a 
subhologram  is  continuous  but  across  subholograra  boundaries  need  not  be 
continuous.  The  format  in  which  the  subholograms  appear  depends  critically 
on  the  style  of  transformation  that  is  being  attempted.  For  example 
Cederqulsc  and  Tal^^^  discuss  the  mapping  of  rings  in  the  (x,y)  plane  to 
point  along  the  v  axis  in  the  (u,v)  plane. 


Clearly  this  transformation  does  not  obey  equation  (4)  in  the  main  text 
and  the  phase  is  not  continuous.  The  hologram  was  therefore  divided  into  N 
concentric  rings  of  average  radius  r^  and  width  Ar  with 


4i(x,y)  »  2xrnx/Xf 


for 

1  Vx^+y^  -  r^  !  <  Ar/2 

A  similar  transformation  can  be  realised  using  the  Mellln  transform  filter 
described  in  the  text  and  a  second  hologram  which  transfers  lines  into 
points.  Both  of  these  hologram  have  continuous  phase  description  and  can 
potentially  have  a  higher  space-bandwidth  product. 

Other  optical  map  transformations  have  been  investigated  by  Bryngdahl^^^ . 
The  first  case  he  considers  is  a  preservation  of  geometric  similarity 

0  •  C)  ■  C) 

Use  of  equations  (3)  result  in  a  phase  function 

<|i(x,y)  =  ~  (x^+y^) 


which  is,  of  course,  the  phase  factor  introduced  by  a  lens  of  focal  length 
f,  ie  a  Fresnel  zone-plate. 


A  second  transformation  described  by  Bryngdahl  used  the  example 

rp  exp(-x/q)  cos(y/q)\ 
kp  exp(-x/q)  sln(y/q)/ 


C)  •  C)  ■  c 


This  transformation  maps  horizontal  lines  Into  circles  of  given  radius 
and  vertical  lines  Into  radii.  The  phase  function  had  the  form 


(J)(x,y)  «  -  exp(-x/q)  cos(y/q) 

The  third  transformation  described  by  Bryngdahl  Involved  using  a  chirped 
one  dimensional  grid  structure  for  local  modification  of  an  Image.  Sections 
of  an  Image  could  be  extended  to  position  relative  to  the  features  In  the 
CGH. 
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APPENDIX  B 


PARAMETERS  AND  PROGRAMS  USED  TO  DEFINE  THE  MELLIN  TRANSFORM  FILTERS 
DESCRIBED  IN  THE  TEXT 


4>(x,y) 


H(-x)  -  1  X  <  0 

H(-x)  -  0  X  >  0 

A  »  514.5  nm 

f  «  1000  mm 

s  *  5/2x  (for  Mellln  Nos  1,3) 

s  =  l/2n  (for  Mellln  No  2) 

4x  *  Ay  =  102.9  ym 

c  *  0.5  (see  Figure  1) 

Anm  =1  (r  <  1  mm) 

Anm  “  (r  >  1  mm) 

r(max)  ■  10.29  mm 

Type  of  CGH  -  type  I  (for  Mellln  Nos  1,2) 

type  II  (for  Mellln  No  3)  (see  Figure  1) 


Bl:  Mellln  No  1 


5  '  nEL'-IN  H  I 
plotter  is  765 
20  P.ftD 

30  LIMIT  0 . 250.0- 130 
40  LOCATE  0- 10e»PRTIO-e.- 100 
50  SCALE  -<12»RATI0>  ..  121:RATI0. - 
12-12 

60  GO SUB  1000 
70  PAUSE 

1000  '  PLOT=  DIM  IM  mm 
1010  Xl=.102?  e  !  X  CELL  SIZE 
1020  Yl=  1029  e  !  Y  CELL  SIZE 
1030  Hl=  0026  e  !  MIN  HEIGHT  OF 
SLIT 

1040  L=  0005145  e  '  WAVELENGTH  0 
F  LIGHT 

1050  F=1000  e  '  FOCAL  LENGHT  OF 
LENS 

•  10^:0  C=  5  e  !  WIDTH  OF  SLIT 
1030  R:S=1  e  !  MIN  RAD  FOR  AMP  CO 

RR 

♦  1075  S=5/'<21:PI>  6  !  SCALE  FACTOR 
FOP  PHASE  CALCULATION 
10S0  FOR  X2=-100  TO  100 
1090  X=X1*X2 

1100  FOR  y2=-IP<SQR'' 10000-X2tX2> 
+1)  TO  IP<SCR<ie0e0-X2*X2)+ 
1) 

3110  Y=yi*Y2 

1120  IF  X2=0  AND  Y2=0  THEN  1170 

1130  GOSUB  2000 

1140  ic^  A=0  THEN  1170 

1150  GOSUB  3000 

1160  GOSUB  4000 

1 170  NEXT  Y2 

! 1S0  NEXT  X2 

1190  RETURN 

2000  '  AMPLITUDE 

2010  P=S0P<X''2+Y'^2) 

2020  IF  R>R9*Y1-'H1  THEN  A=0  0  GO 
TO  2050 

2030  IF  P>R9  THEN  AsRS^R  8  GOTO 
2050 

2040  A=1 
2050  RETURN 

7000  !  PHASECREF  SAITOH 
3010  !•=■  X2=e  THEN  P=ABS';Y)CPI^2 
e  GOTO  3040 

3020  P=Y»ATN<Y/X>-  .5»X»L0G<X''2  +  Y 
^2>+X 

3030  IF  X<0  THEN  P=P+Y»PI 
3040  P=-<:2»PI»P»6.''<L*F>  > 

3050  P=ATN2'!'SIN<P) -COS<P)  ;■ 

3060  RETURN 

♦  4000  •  HOLBAR 
/v4ei0  Y5=A»Y1».5 
^4020  C1=X+P*X1-'<2»PI> 

04030  Y8=Y+Y5 

0  4040  V9=Y-Y5 
^4060  X8»Cl-XI-"4 


.  4070 
04110 
♦  4120 
5000 
5010 
5020 
5030 
5040 
505© 

5051 

5052 

5053 

5054 

5055 
5060 


X9=C1+X1 M 
GOSUB  5000 
RETURN 
MOVE  X8-Y9 
DRAW  X9-Y9' 

DRAW  X9-Y8 
DRAW  XB-YB 
DRAW  X8-Y9 
PENUP 

FOR  X7=X8  TO  X9  STEP 

MOVE  X7-Y9 

DRAW  X7-Y8 

NEXT  X7 

PENUP 

RETURN 


.01 
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B2:  Mellln  No  2 


5  '  MELLIN  #  2 

4070 

10  PLOTTER  IS  70 

5 

V  41  10 

20  RAD 

♦  4120 

30  LIMIT  0.250.0 

,  180 

5000 

40  LOCATE  0,  1004RATICI.0..  100 

sni  A 

50  SCALE  -<:i21:RATI0::>,  12»RRT10^- 

5020 

17.12 

5030 

^ci  cnsuB  I00fl 

5040 

70  PAUSE 

5050 

1000  '  PLOT:  dim 

I  N  mm 

5051 

1010  Xl= . 1029  e  ' 

X  CELL  SIZE 

5052 

1020  Yl=. 1029  e  ! 

Y  CELL  SIZE 

5053 

1030  HI =0026  e  ! 

MIN  HEIGHT  OF 

5054 

SLIT 

5055 

1040  L  =  .  0005145  e  .'  WAVELENGTH  0  5068 

F  light 

1050  F=1000  e  !  FOCAL  LENGHT  OF 
LEM^ 

10t‘0  C.=  .5  e  ‘  WIDTH  OF  SLIT 
ie"0  R9=l  §  !  MIN  RAD  FOR  AMR  CO 
FR 

♦  1075  S=1/''2*PI>  @  !  SCALE  FACTOR 
FOP  PHASE  CALCULATION 
10S0  FOR  X2=-100  TO  100 
1090  X=X1*K2 

1100  FOR  Y2=-IP<SQR<  18000-X2*.v:2) 

+1>  TO  IP(SQR<10e0e-X2»X2>+ 

1  > 

1110  Y=Y1»Y2 

1120  IF  X2=0  AND  V2=0  THEN  1170 

1130  GOSUB  2000 

1140  IF  A=0  THEN  1170 

1150  GOSUB  3000 

1160  GOSUB  4000 

1170  NEXT  Y2 

11S0  NEXT  X2 

1190  RETURN 

2000  !  AMPLITUDE 

2010  P=S0R<X'^2+Y'^2) 

2029  IF  R:  R9»Y1''H1  THEN  A=0  0  GO 
TO  2050 

2030  IF  PR?  THEN  A=R9/R  0  GOTO 
2050 

2040  A=1 
2050  RETURN 

3000  '  PHASE ‘  REF  SR I  TO > 

3010  IF  X2=0  THEN  P=ABS<y )*PIx2 
e  GOTO  3040 

3020  P=Y»ATN< V  'X)  -  .St^tLOGcX'-'E-^Y 
''2)  +  :< 

3030  IF  X^0  THEN  P=P+Y»PI 
3040  p=-<'2*PI»:P»£/'<L»F)  .> 

3050  P=ATN2<SIN<P) ,COS<P) ) 

3060  RETURN 
♦4000  !  HOLBAR 
^4010  Y5=flTYlf.5 
14020  C1  =  X+P»X1 '(2»PI> 

V403e  VS=Y-fY5 
04040  V?*Y-V5 
^4060  X8=Cl-Xl/4 


X?=C ! +X1  ^  4 
GOSUB  5000 
RETURN 
MOV/E  X8>Y9 
DRAW  X?,Y? 

DRAW  X9,Y8 
DRAW  X8,Y8 
DRAW  X8,Y9 
RENUP 

FOR  X7=X9  TO  X9  STEP 

MOVE  X7.Y9 

DRAW  X7,Y8 

NEXT  X7 

PENUP 

RETURN 


.  01 
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B3:  Mellin  No  3 


5  '  MELLIN  «  3 
le  PLOTTER  IS  765 
26  RftD 

36  LIMIT  0.250.0, 186 
40  LOCPTE  e,  ie6»RflT  10,0 . 100 
50  SCftLE  -< 12«RftTI0>, 12»RflTI0,- 

4  ^  ^ 

^  ^  f  X  w 

6C’  OGiUB  1000 
70  PftUSE 

1000  •  PLOT:  DIM  IN  mm 
1010  Xl=.ie29  e  X  CELL  SIZE 
1020  Yl=  1029  @  I  Y  CELL  SIZE 
1036  hl=.00£6  6  <  MIN  HEIGHT  OF 
SLIT 

1040  L=. 0005145  e  ?  WAVELENGTH  0 
F  LIGHT 

1050  F=1000  e  !  FOCAL  LENGHT  OF 
LENS 

1000  C=  5  e  !  WIDTH  OF  SLIT 
1070  R9=l  0  !  MIN  RAD  FOR  AMP  CO 
RR 

♦  1075  S=5,^'''2*PI>  0  !  SCALE  FACTOR 
FOP  PHASE  CALCULATION 
1080  FOR  X2=-100  TO  100 
1090  X=XlxX2 

1100  FOR  Y2=-IPCSQR<  lOeOO-XZCXZ') 
♦1>  TO  IP<SQR(10000-X2»X2l+ 
1> 

111©  V=V1»V2 

1120  IF  X2=0  AND  Y2=0  THEN  1176 

1130  GOSUB  2600 

1140  IF  A=0  THEN  1170 

1150  GOSUB  3000 

1160  GOSUB  4000 

1176  NEXT  Y2 

1180  NEXT  X2 

1190  RETURN 

2000  !  AMPLITUDE 

2010  R=S0R<X'*2+Y''2) 

2020  IF  P>R9*YI/'H1  THEN  A=0  0  GO 
TO  2056 

2030  IP  R' R9  THEN  A^RS/R  0  GOTO 
2050 

2040  A=1 
2050  RETURN 

3000  «  PHASECREF  SAITO) 

3010  IF  X2=0  THEN  P=ABS<; Y)*PI/'2 
6  GOTO  3040 

3020  P=y»ATN <  V/X > - . 5*X»L0G  < X ' 2  +  Y 
'•2)+X 

3030  IF  X<0  THEN  P=P+Y»Pr 
3040  P=-<2*PI»P»S/'(L»F>) 

3050  P=ATN2<:SIN<:P),C0S<P>) 

3060  RETURN 
#4000  !  H0L8AR 
^4010  Y5=A4Y1*.5 
X4020  C1=X+P*X1/-<2»PI> 

04030  V9=V+V5 
04040  Y9»Y-Y5 

^4050  IF  P<-<PI/2>  THEN  GOTO  4680 


A  4060 
^4070 

0  4080 
0 

i 

0 

0^109 

0. 

0 

4^120 

5000 

5005 

5010 

5020 

5030 

5040 

5050 

5051 


4090 


4110 


5052 

5053 

5054 

5055 
5060 


IF  P-PI.2  THEN  GOTO  4166 
X8=C1-X1  '4  0  X9=C1+X1.'4  0  G 
OSUB  5000  0  RETURN 
X8=X-Xl/'2  0  X9=Cl+Xl/4  0  CO 
SUB  5000 

X8=Cl+3*Xl/4  0  X9=X+Xl/2  0 
GOSUB  5060  0  RETURN 
X3=C1-X1.'4  0  X9=X+XI/'2  0  GO 
SUB  5000 

X3=X-X1^'2  0  X9=Cl-3»Xl/4  0 
GOSUB  5006  0  RETURN 
RETURN 

IP  X9-X8<H1  THEN  5066 

MOVE  Xe,Y9 

DRAW  X9,Y9 

DRAW  X9,Ye 

DRAW  X8,Y8 

DRAW  Xe,Y9 

RENUP 

FOR  X7=XS  TO  X9  STEP  .01 

MOVE  X7,Y9 

DRAW  X7,y8 

NEXT  X7 

PENUP 

RETURN 
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Abstract 

The  procedure  for  designing  a  detour  phase  computer  generated  holograms  is 
described  with  particular  reference  to  a  scale  and  rotation  invariant  transforma¬ 
tion  filter  known  as  a  Mellin  transform  filter.  Experimental  results  are 
presented  which  demonstrate  the  main  features  of  the  devices. 


